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ABSTRACT 

We have carried out a study of the chemical abundances of ^H, '^He, ^^C, ^^C, ^''N, ^^N, ^^O, 
i7q^ 20]\jg g^jjjj ^^p^g symbiotic stars (SSs) by means of a population synthesis code. We find 
that the ratios of the number of 0-rich SSs to that of C-rich SSs in our simulations are between 
3.4 and 24.1, depending on the third dredge-up efficiency A and the terminal velocity of the stellar 
wind v{oo). The fraction of SSs with extrinsic C-rich cool giants in C-rich cool giants ranges 
from 2.1% to 22.7%, depending on A, the common envelope algorithm and the mass-loss rate. 
Compared with the observations, the distributions of the relative abundances of ^^C/^'^C vs. 
[C/H] of the cool giants in SSs suggest that the thermohaline mixing in low-mass stars may exist. 
The distributions of the relative abundances of C/N vs. 0/N, Ne/0 vs. N/0 and Hc/H vs. N/0 
in the symbiotic nebulae indicate that it is quite common that the nebular chemical abundances 
in SSs are modified by the ejected materials from the hot components. Helium overabundance 
in some symbiotic nebulae may be relevant to a helium layer on the surfaces of white dwarf 
accretors. 

Subject headings: binaries: symbiotic — accretion — stars: AGB — Galaxy: stellar content 



1. Introduction 

Symbiotic stars (SSs) are usually interacting bi- 
naries, composed of a cool star, a hot component 
and a nebula. The cool component is a red giant 
which is a first giant branch (FGB) or an asymp- 
totic giant branch (AGB) star. The hot compo- 
nent is a white dwarf (WD), a subdwarf, an accret- 
ing low-mass main-se q uence s tar, or a neutr o n star 
(iKenvon fc WebbinkI [iQSl iMiirset et all Il991 



quiescent phase or the outburst phase. During 
the quiescent phase, SSs are undergoing a stable 
hydrogen burning on the surfaces of the WD ac- 
cretors. The outbursts may be due to either the 
thermonuclear outbursts on the surfaces of ac- 
creting WDs which are called as symbiotic novae 
( Tutukov fc Yungelson 19761 : Paczvhski fc Rudak 
1980f). or accretion-disk instabilities around a 
nearly steady burning WDs which are called 



Yungelson et al.lll995t1lben fc Tutukovlll996l ). In 

general, SSs are low-mass binaries with an evolved 
giant transferring materials to a hot WD compan- 
ion which bu rns the accreted m aterials more or 
less steadily (jMikolaiewskal 120071 ) . According to 
the observed characteristics of their specular spec- 
tra and photometries, SSs may stay in either the 
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as multiple outbursts ijMitsumoto et al 
^A-Iikolajcwska 2007,)- Recent re views of the prop- 



erties of SSs can be found in Miirset fc Schmid 
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(Il999l) and lMikolaiewskal (I200ll2007l) . 

The cool components usually have a high mass- 
loss rate in the SSs with WD accretors. Their 
stellar winds have the chemical abundances of the 
red giants. Parts of the stellar winds are ionized 
by the hot components. The abundance deter- 
minations with nebular diagnostic tools are possi- 
ble (S child. Boyle fc Schmid ..1992 ). SSs provide a 
good opportunity for measuring the abundances of 
the red giants. The hot components may support 
an additional high-velocity wind during the sym- 
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biotic nova outbursts (jKenvon fc Webbinklll984[ ). 
High-velocity outflows have been observed in very 
broad e missio n lines i n essentia lly all the symbiotic 
novae (jKennv fc TavloJ l2005l ). iLivio fc TruTM 
(jl994l ) analyzed classical nova abundances and 
concluded that enhanced concentrations of heavy 
elements are significant. Therefore, the abun- 
dances of the ejected materials from the hot com- 
ponents in SSs may be different from those of the 
cool giant stellar winds. However they are simi- 
lar to those of the classical novae. In general, the 
hot components have high luminosity and effective 
temperature (They are similar to th e central stars 
of pla netary nebulae. See Figure 5 in lMiirset et al.l 
so that they can ionize the nebulae. In 
some eruptive SSs, the nebulae can also be ion- 
ized by the region where the winds from t he hot 
and cold components collide (jWillson et al ]|i984). 
It is difficult to determine the origin of the nebu- 
lae in SSs. Based on emi ssion line fluxes from Cm , 
Civ, Niii, Njv and OTTT. lNussbaumer et al. I (|l988l) 
found that the CNO abundance ratios of the neb- 
ulae in SSs are in the transition region from gi- 
ants to supergiants and concluded that the nebula 
is mainly due t o mass lost by the red giant . On 
the other hand, iNussbaumer fc Vogel (1 19891) sug- 
gested that the flux ratios of the emission line de- 
pend on the relative abunda nces of the two winds 
i n the ir study on Z And. IVogel fc Nussbaumer 



(11992) found that CNO abundance ratios of PU 
Vul are of the characteristic of novae. 

In short, SSs offer an exciting laboratory for 
studying novae, the red giants and the interaction 
of the two winds. The chemical abundances of 
the ejected materials from the hot components, 
the stellar wind of the cool components and the 
symbiotic nebulae are the key factor to under- 
stand them. Up to now, a series of observa- 
tional data of the chemical abundances in SSs 
have been published 



ONe WD accretors. This makes it possible to 
simulate the chemical abundances of the sym- 
biotic novae. [Grocncweecn fc de Jong (19 9^), 



Wagenhuber fc GroenewegenI (Il99i_.. 
(I2OO2I). Ilzzard et al.' h004). Ilzzardi (|2004l ) and 
Marigo fc Girardi (.2007.) developed thermally 
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pulsing asymptotic giant branch (TP-AGB) syn- 
thesis by which the chemical abu ndances of the 

cool co mponents can be computed. iLii. Yungelson fc Hani 
( 20061 ) (hereafter Paper I) constructed the models 
for the SSs' populations. According to the above 
results or models, it is possible to set up a prelim- 
inary model for theoretical study on the chemical 
abundances of SSs. 

In the present paper, the study of the chemical 
abundances of SSs are carried out by means of 
a population synthesis code. In § 2 we present 
our assumptions and describe some details of the 
modeling algorithm. In § 3 we discuss the main 
results and the effects of different parameters. § 4 
contains the main conclusions. 

2. Model 

For the binary evolution, we use a rapid binar y 
star evolution (BSE) code of iHurlev etal] (|2002t ). 
Below we describe our algorithm from several as- 
pects. 

2.1. Symbiotic Stars 

Paper I carried out a detailed study of SSs. 
Comparing the observed distributions with the 
those predicted of the orbital periods and the hot 
component masses in SSs, it is worth to mention 
that there is remarkable disagreement between the 
observations and predictions. 



Mikolaiewskal (|2007t ) showed the distributions 



Nussbaumer et al. (1988) 



De Freitas Pa checo fc CostaVl992'). ICosta fc de Freitas 
( 1994 ) andlLuna fc Costa (2005) gave th e chemi- 



of the measured orbital period in about 70 SSs. 
The vast majority (~ 75%) of these systems have 
— orbita J periods shorter than 1000 days while Paper 
cted the distribution of longer orbital peri- 
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cal ab undanc es of some symb i otic ne bulae. ISchild. Bovle fc Schmid 
(|l992l ) and ISchmidt et aP (|2006l ) showed the 
chemical abundances of the cool giants in sev- 
eral SSs. However, there are few theoretical 
studies about it 



odB. Tha main reasons results from two aspects: 



It is hard to measure long orbital periods 
A ll the SSs with measur ed o rbital periods 



. . Recently , iKovetz fc Prialnik 

(1997) and Jose fc Hernanz ( 19981 ) carried out 
detailed numerical simulations for the chemi- 
cal abundances of the novae with CO WD or 



— .Belczvnski et al. I (l2000l ) and I Mikola i e wskal 
(I2OO3I) are S-type and about half of them are 
eclipsing binaries. Considering the ampli- 
tude of radial velocity changes or eclipses, 
it is easier to measure the short orbital 



2 



periods. Ab o ut 15 percent of 27 SSs in 
Mikolajewskal ( 2003 ) have longer orbital pe- 



riods th an 1000 days while it is about 25 per- 
cent in iMikolaiewskal (|2007l ) . There would 
be more SSs with long orbital periods with 
further detailed observations. 

2. In Paper I, SSs are detached binary sys- 
tems and the process of mass transfer re- 
sults from the stellar wind of cool giants. 
To our knowledge, if the mass ratio of the 
components {q = Mdonor/-^accretor) at on- 
set of Roche lobe overflow (RLOF) is larger 
than a certain critical value Qc, the mass 
transfer is dynamically unstable and results 
in the formation of a common envelop in 
short time scale (about 100 years). The is- 
sue of the criterion for d ynamically unsta- 
ble R LOF Qc is stiU open. iHan et al.l (|2001 
20021 ) showed that Qc depends heavily on 
the assumed mass-tra nsfer efRciency . Pa - 
per I adopted qc after Hurley et al. ( 2002h 



in which the mass-transfer efficiency is 1. In 
Paper I, q^ is usually smaller than 1.0. The 
observed MgiantA^wos in SSs are between 
2.0 and 4.0 and the predicted those in Paper 
I are from 1.0 to 4.0, which means that there 
should not be stable RLOF in SSs. However, 
Mikolaiewskal (|l007) suggested that RLOF 



can be quite common in symbiotic binaries 
with orbital periods shorter than 1000 days. 
If this were true, the theoretical model of 
mass transfer should be advanced. But this 
is out of the scope of this paper. 

The disagreement between the observed and the 
predicted distributions of orbital periods in Pa- 
per I resulted from observational biases and poor 
knowledge of the mass transfer mechanism. 



Mikolaiewskal (|2007f) showed that most of SSs 
have WD masses less than 0.6 Mq, while Paper 
I predicted the distribution of higher WD masses. 
The disagreement is due to the following: 

1. For mass estimates, an orbital inclination 
oi i = 90" o r a lim it to i (see table 2 in 
iMikolajewskal (|2003h ) is typically assumed, 
hence the estimates are lower limits. 

2. Pap er I use d the ma ss-loss law suggested by 
IVassiliad is fc WoodI (.199 3) for AGB stars. 
However, Mikolaiewskal ( 20071 ) suggested 



that both the symbiotic giants and Miras 
have higher mass-loss rates than single gi- 
ants or field Miras, respectively. Therefore, 
Paper I may have overestimated the hot 
component masses. 

In short, due to poor knowledge of the mass 
loss from the giants and the mass transfer mecha- 
nism in SSs, Paper I can only crudely predict some 
observed characteristics. In this work, we accept 
still all the criterions and the concepts on SSs in it. 
Following Paper I, we assume that binary systems 
are considered as SSs if they satisfy the following 
conditions: 

• The systems are detached. 



The luminosity of the hot component is 
higher than lOL© which is the 'thresh- 
old' luminosity fo r the hot conipone nt of 
SSs as inferred by iMiirset et al.l (|l99ll ) and 
Mikolaiewska fc Kenvon ( 1992f ). This may 
be due to the thermonuclear burning (in- 
cluding novae outbursts, stationary burning 
and post-eruption burning). 



• The hot component is a WD and the cool 
component is a FGB or an AGB star. 

The liberation of gravitational energy by the ac- 
creted matter may make the luminosity of the 
component larger than lOL©. A detailed accre- 
tion model of SSs was discussed in Paper I. Here, 
we do not model it. 

All symbiotic phenomena in our work are pro- 
duced by the hydrogen burning on the WD sur- 
face. All SSs should stay in one of the following 
three phases: 

• the stable hydrogen burning phase; 

• the thermonuclear outburst phase (symbi- 
otic nova); 

• the declining phase after a thermonuclear 
outburst. 

Paper I and the present paper do not simulate 
the accretion disk in SSs. As mentioned in the 
Introduction, the multiple outbursts due to the 
accretion-disk instabilities u sually occur aroun d 
nearly steady burning WDs (jMikolaiewskall2007t ). 
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In this work, they can be included in the quies- 
cent SSs which are undergoing the stable hydrogen 
burning. 

SSs are complex binary systems. There are 
many uncertain physical parameters which can 
affect the population of SSs. Paper I showed that 
the numbers of SSs and the occurrences of symbi- 
otic novae are greatly affected by the algorithm of 
common envelope evolution, the terminal veloc- 
ity of stellar wind v{oo) and the critical ignition 
mass AM^^ which is necessary mass accreted by 
WDs for a thermonuclear runaway. The struc- 
ture factor of the stellar wind velocity aw and an 
optically thick wind give a small uncertainty. In 
this work, we use the models of SSs in Paper I 
and discuss the effects of the algorithm of com- 
mon envelope evolution, v{oo) and AM^^P on the 
chemical abundances of SSs. Other subordinate 
parameters are the same as in the standard model 
of Paper I. 

Common Envelope: For the common envelope 
evolution, it is generally assumed that the orbital 
energy of the binary is used to expel the envelope 
of the donor with an efficiency acc- 



bind 



orb ; 



(1) 



where i?bind is the total binding energy of the en- 
velope and A Eorh is the orbit a l ener gy released in 
the spiral-in. Nelemans et al. ( 2000l ) suggested to 
describe the variation of the separation of compo- 
nents in the common envelopes by an algorithm. 
This algorithm is founded on the equation for the 
system orbital angular momentum balance which 
implicitly assumes the conservation of energy: 



AJ _ Mo 



(2) 



where J is the total angular momentum and A J is 
the change of the total angular momentum during 
common envelope phase. In the above formula, M 
and M(. are respectively the mass and the envelope 
mass of t he donor, and m is t he com panion mass. 
Following iNelemans fc Tout! (|2005h and Paper I, 
we call the formalism of Eq. ^ a-algorithm and 
that of Eq. ^ 7-algorithm, which are respec- 
tively simulated in different cases (see Tabic [T|) . 
We take the 'combined' parameter accAco as 0.5 
for a-algorithm. Acc is a structure parameter that 
depends on the evolutionary stage of the donor. 



For 7-algorithm, 7 = 1.75. 

v{oo): It is difficult to determi ne the term! 
nal ve locity of stellar wind ^(oo). IWinters et al 



fitted the relation between the mass-loss 
rates and the terminal wind velocities derived from 
their CO (2-1) observation by 

logio(Af/Af0yr-i) = -7.40+^ logio(z;(oo)/kms-i). 

(3) 

T he mass-loss rate is given by the formulation 
of IVassiliadis fc WoodI (jl993h or iBlockeil (|l995l) 
and v{oo) can be obtained by Eq. ([3]). How- 
ever, Eq. ([3]) is vahd for M close to IQ-^M© 
yr~^. For a mass-loss rate higher than 10~^ Mq 
yr~^, Eq. Q gives too hig h v(oo ). Based on 
the models of Winters et al.l ( 200Cl[ ). we assume 
v{oo) — 30kms~^ if v{oo) > BOkms""'^. In the 
standard model of Paper I, v{oo) = ^Wcsc where 
Vcsc is the escape velocity. In the present paper, 
we also carry out various v{oo)s simulations. 



AMj^P: The critical ignition mass of the no- 
vae depends mainly on the mass of accreting 
WD, its tempera t ure a nd material accreting rate. 
Yungelson et al. 1 (|l995h gave the 'constant pres- 
sure' expression for AM^f* as 



AMWP 

^ = 2 X 10" 

Ma, 



f AL 



WD 



WD 



-0.8 



(4) 



where Mwd is the mass of WD accretor and i?wD 
is the radius of z ero-temperature degenerate ob- 
jects (INauenberg 1972), 



i?WD = O.O112i?0[(MwD/Mch 



,-2/3_ 



where M, 



ch 



■(MwD/Mch) 
(5) 

I.433M0 and i?0 = 7 X 10^° cm. 



2/31I/2 



1,0 



Nelson et al.l (|2004l ) gave numerical fits to the crit- 
ical ignition masses f or nov ae models calculated by 



Prialnik fc Kovetz 



( T995I) . In most simulations 
we adopt Eq. ^ for AM^°. However, in or- 
der to investigate the influences of the AM^-^ on 
our res ults, we also c a rry ou t a simulation for Eq. 



(Al) of iNelson et al.l (|2004[ ) in which AM<^° de- 



pends on the mass accretion rates and masses of 
WD accretors. 
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2.2. Chemical Abundances on the Surface 
of the Giant Stars 

For a single star, three dredge-up processes 
and hot bottom burning in a star with initial 
mass higher than 4Mq may change the chemical 
abundances of the stellar surface. We accept the 
prescr iptions of §3.1.2 and 3.1.3 in llzzard et al 



(|2004l ) for the first dredge-up during the first 
giant branch and the second dredge-up during 
early asymptotic giant branch (E-AGB). For the 
third dredge-up (TDU) and the hot bottom burn- 
ing during the TP-AG B phas e, we use the TP - 
AGB synth esis inlGroen ewcgc n fc de Jone j 1993 
iKarakas. Lattanzio &: Fold (|2002|). Izzard et al 



()2004l ). llzzardl()2004f ) and lMarieo fc Girardil()2007l) . 
All details can be found in Appendix A. In the 
present paper, we give the chemical evolutions of 
iH, 4He, 12c, "C, "N, i^N, i^o, i^q, 20^^ 
22 Ne on the stellar surface. 

SSs in our work are binary systems. The binary 
mass transfer can change the chemical abundances 
of the stellar surface. In binary systems, there are 
two ways to transfer mass: (i)accretion from the 
stellar wind material of a companion; (ii)Roche 
lobe overflow. BSE model cont ains a standard 
Bond i-Hoyle type wind accretion ( Bondi fc Hovld 



19441 ) and a conservative mass transfer during the 
stable Roche lobe overflow. After a star obtains 
AM from its companion, the chemical abundances 
of the stellar surface X2 is given by 



void ^ /\/fcnv 



Al X AM 



AM 



(6) 



where M°"'^ is the envelope mass and Xi is the 
chemical abundances of its companion. 

2.3. Chemical Abundances of the Ejected 
Materials from the Hot WD Accre- 
tors 

The hot WD accretors eject ma t erials from 



their su rfaces. iKovetz fc PrialnikI (|l997l ) and 



[jose fc Hernanj ( 1998h showed detailed study 



their chemical abundances during the thermonu- 
clear outbursts. There is no direct observational 
evidence that SSs in the stable hydrogen burn- 
ing phase and the declining phase lead to winds 
from the hot components. Theoretically, the hot 
companions have high luminosity and tempera- 
ture during the above two phases. Some parts 



of the materials on the surface of the accreting 
WDs may be blown o fi' due to high luminosity 
(llben fc Tutukovlll99(ih . To our knowledge, no 
literature refers to them in detail yet. However, 
their chemical abundances are possibly between 
those of the ejected materials during the ther- 
monuclear outbursts and the stellar winds from 
the giant stars. In the present paper, we only give 
the chemical abundances of the ejected materials 
during the thermonuclear outbursts. 

Based on the typical characteristics of the ob- 
served systems, nova outbursts are usually di- 
vided into symbiotic novae and classical novae. 
The classical novae can only last several days or 
weeks and the variation of their visual magni- 
tudes are between -5 and -10. The symbiotic 
novae usually last several decades and the vari- 
ation of their visual magnitudes are between -3 
and -8. Thermonuclear runaways appear to be 
the most promisin g mechanis m for classi c al no- 
vae (|Kenvonlll986h . [Tutu kov fc YungelsonI (|l976f ) 
suggested that the above mechanism is applicable 
to symbiotic novae. The observational differences 
between the symbiotic novae and the classical no- 
vae may stem from the prope rties of the binary 
systems ( Iben fc Tutukovl 1996f) . The physical na- 
ture of the classical novae and the symbiotic no- 
vae is identical. Therefore, we assume that the 
thermonuclear outbursts occurring in our binary 
systems selected for SSs are the symbiotic novae. 

Numerical studies of the nova outbursts have 
been carried out by Shara. Prialnik fc Kovetz 



1993 
1995 



Covetz fc PrialnikI ( 19941) :lPrialnik fc Kovetz 
Yaron et all (|2005L iKovetz fc PrialnikI 



19971 ) published detailed multicycle calculations 



of the nova outbursts for CO WDs with mass 
ranging from 0.65 to 1.4 Mq. In their simula- 
tions, the element abundances of nova ejecta are 
affected by the four basic and independent pa- 
rameters: C/0 ratio in the accreting WD, its 
core temperature Twd, the mass accreting rate 
MwD and the mass Mwd- In order to test the ef- 
fect of the WD composition on t he ab undances of 
nova ejecta, iKovetz fc PrialnikI (|l997l ) calculated 
the models with the accreting WDs composed of 
pure-C, pure-0 and C/0=1, respectively. They 
found that the WD composition is not reflected in 
the abundances of ejecta. In an extensive study 
of close binary evolution, llben fc Tutukovl ( 19851 ) 



showed that the AGB phase may be suppressed in 
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a close binary and a WD formed in such a system 
should have a ratio very close to unity. In our 
work, we assume that the CO WD composition is 
C/0=1. A WD temp erature of 10^ K corresp onds 
to an age of 10^ yr (jlben k Tutukovl 119841) . In 
fact, for most of SSs, the time interval between 
the formation of the WD and the beginning of 
the symbiotic sta ge may be lon ger than it, up to 



10^" yr (Yungc lson et al.Ml995r ). Neglecting the 



effect of the nova outbursts on the temperature 
of WD accretor Two, we assume that it is 10'' K 
We select 10 models in iKovetz fc Prialnid (|l997h 
in which Two = 10'' K and C/0 ratio of the ac- 
creting WD is 1. Their element abundances are 
determined by the mass accreting ra te Mwd and 
Mwp. By a bilinear interpolation dPress et al. " 



19921 ) of 10 models m 
the abundances of ^H. 



Kovetz fc PrialnikI (|l997n . 

4He, 12c, 13c, "N, 15N, 



1^0 and in the nova ejecta are calculated. If 
Mwd or Mwd in SSs are not in the range of the 
bilinear interpolation, they are taken as the most 
vicinal those in the 10 models. 

Kovetz fc PrialnikI (jl997l ) di d not give the 
abundances of ^^Ne and ^^Ne. Jose fc Hernand 



(|1998[ ) gave the nucleosynthesis in the nova out- 
bursts with CO and ONe WDs. In their calcu- 
lations, the nova nucleosynthesis are affected by 
Mwd, Mwd, the initial luminosity (or Twd) and 
the degree of mixing between core and envelope. 
To test the effect of the WD mass, they carried 
out a number of simulations involving both CO 
WD {Mwd=0.8, 1.0 and 1.15 Mq) and ONe ones 
(MwD=l-0, 1.15, 1.25 and 1.35Mo). Mwd is 
2 X lO^i^Af© yr~i and their initial luminosity is 
\Q~^Lq. The degree of mixing between the core 
and the envelop e is a very uncertain parameter. 
Jose fc Hernand ( 19981) modeled three different 
mixing levels: 25%, 50% and 75%. Their results 
showed that higher mixing degree favors the syn- 
the sis of higher metal nuc lei in ONe WDs. Follow- 
ing Starrfield et al. ( 1998h . we adopt a 50% degree 
of the mixing. F o r the simulations of CO WDs in 
Jose fc Hernand (1998), we select the three nova 
models (The degree of the mixing is 50%) to cal- 
culate the abundances of ^''Ne and ^^Ne by the 
fitting formulae: 



logX(20Ne) = -3.206 + 0.14476A-%d, 
logX(22Ne) = -2.57118 + 0.60784Mwd 

(7) 

where Mwd is in solar unit and the formulae agree 



with the numerical results to within a factor of 1.1. 
In the above fits, the abundances of ^^Ne and ^^Ne 
depend weakly on Mwd- The range of CO WD 
mass Mwd is from ^ 0.5 to 1.4 M©. The abun- 
dances of ^''Ne calculated by Eq. ([7]) are between 
~ 10-3-13 and 10-3 00, and the abundances of ^^^e 
are between ~ iQ-^-SS ^^^^ 10-3-38_ Therefore, we 
use Eq. to calculate X(20Ne) and ^(^^Ne) 
of all symbiotic novae with the CO WDs. For 
the nova of ON e WD, we fit data of Table 3 in 
Jose fc Hernand (1998) for a 50% degree of mix- 
ing by: 



-1.486 + 1.982MwD - 0.992M^d, 
-0.839 - 0.103MwD + 0.197M2^j3' 



-10.664+ 14.798A/WD 
-14.691 + 22. 517A/WD 
4.504- 11. 174MwD 
-2.196 - 2.313MwD + 2.404M^j3, 



6.025M2,j3, 
9.607M2,j3, 
5.079M2,j3, 



-9.835 + 17.609Mwd 
-14.492 
-0.536 - 0.192MwD 



•550M2,j„ 
22.095MwD - 9.366M^o, 



-21.20356+ 34.46454MwD - 15.9763M^d, 



the numerical results to within 



0.33769M2^ 



logX(iH) = 
logX(4He) = 
logX(i2c) = 
logX(i3c) = 
logX(i4N) = 
logXpN) = 
logX(ioo) = 
logX(i^O) = 
logX(20Ne) = 
logX(22Ne) = 

which agree with 
a factor of 1.3. 

We neglect other elements because their abun- 
dances are much smaller than the above elements 
or their isotopes. All abundances in the nova 
ejecta are renormalized so that their sum is 1.0. 

2.4. Symbiotic Nebulae 

In SSs, there may be two winds. One comes 
from the cool giant and its chemical abundances 
are similar to those of the red giant envelope and 
its velocity is between ~ 5 and 30 km . The 
other may come from the hot components. The 
components during the symbiotic outbursts have 
the winds with a high velocity 1000 km s-^) 
although there is not a detailed description on the 
wind from the hot components du ring the quies - 
cent phase and the declining phase ( Kenvon|[l986f ) . 
Wind collision is i nevitable in the s y mbiot i c out- 
bursts. Recently, iKennv fc Tavloil (|2005l l2007l ) 
carried out a detailed study for the colliding winds 
in SSs. According to their simulations, the struc- 
ture of the colliding winds (including temperature 
and density) is very complicated. Symbiotic nebu- 
'?ae should have a similar structure with the collid- 
ing winds at least during the symbiotic outbursts. 
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In order to avoid the difficulty of modeling 
the real symbiotic nebulae and study the nebu- 
lar chemical abundances by population synthesis 
method, we assume roughly that the compositions 
of the symbiotic nebulae undergo two independent 
phases: 

• The symbiotic nebula is mainly composed 
of the ejected materials from the hot WD 
when the thermonuclear runaway occurs and 
this phase last for ton- The lasting time 
scale ton on which the hot components are 
in a 'plateau' phase with high luminosity is 
given by Eq. (30) of Paper I. At this phase, 
the symbiotic nebula embodies the chemical 
characteristics of a nova. 

• After ton, the symbiotic nebula is mainly 
composed of the stellar wind materials from 
the cool giants until the next symbiotic nova 
occurs. During this phase, the symbiotic 
nebula embodies the chemical characteristics 
of the cool giants in our models. For stable 
hydrogen burning SSs, the symbiotic nebula 
is always composed of the stellar wind ma- 
terials from the cool giants. 

However, some symbiotic novae have helium 
WD accretors. To our knowledge, no literature 
refers to their chemical abundance. According to 
Paper I, the occurrence rate of the symbiotic no- 
vae with helium WDs is at most about ^ of that 

75 

with CO and ONe WDs. We do not consider SSs 
with helim WD accretors in this paper. 

2.5. Basic Parameters of the Monte Carlo 
Simulation 

We carry out binary population synthesis via 
Monte Carlo simulation technique in order to ob- 
tain the properties of SSs' population. For the 
population synthesis of binary stars, the main in- 
put model parameters are: (i) the initial mass 
function (IMF) of the primaries; (ii) the mass-ratio 
distribution of the binaries; (iii) the distribution of 
orbital separations; (iv) the eccentricity distribu- 
tion; (v) the metallicity Z of the binary systems. 

We use a simp l e app roximation to the IMF 
of iMiller fc Scalol (Il979h . The primary mass 



Eggleton etahl (|l989l ) 



mpl e app ro 
13 (|l979h . 

is generated using the formula suggested by 



Ml 



0.19X 



(l-X)0-75 + 0.032(1 -X)0-25' 



(9) 



where X is a random variable uniformly dis- 
tributed in the range [0,1], and Mi is the primary 
mass from O.SM© to 8Mq. 

For the mass-ratio distribution of binary sys- 
tems, we consider only a constant distribu tion 
(|Mazeh et al.lll992l : IColdberg fc Mazehlll994 . 



n{q)^l, 0<g<l, (10) 

where q = M2/M1. 

The distribution of separations is given by 



loga = 5X + 1, 



(11) 



where X is a random variable uniformly dis- 
tributed in the range [0,1] and a is in Rq. 

In our work, the metallicity Z=0.02 is adopted. 
We assume that all binaries have initially circu- 
lar orbits, and we follow the evolution of both 
components by BSE cod e, including the effe ct of 
tides on binary evolution (jHurley et al.ll2002l ) . We 
take 2 x 10^ initial binary systems for each simu- 
lation. Since we present, for every simulation, the 
results of one run of the code, the numbers given 
are subject to Poisson noise. For simulations with 
2 X 10^ binaries, the relative errors of the numbers 
of the symbiotic systems in different simulations 
are lower than 5%. Thus, 2 x 10^ initial binaries 
appear to be an acceptable sample for our study. 

We assume that one binary with Mi > O.SMq 
is formed annua lly in the Galaxy to calculate the 
birthrate of SSs ( Yungelson et al.l 19941 : Han et al 
1995al lbh. 



3. Results 

We construct a set of models in which we vary 
different input parameters relevant to the chem- 
ical abundances of SSs. Table [T] lists all cases 
considered in this work. Many observational evi- 
dences showed that the terminal velocity of stellar 
wind v (oo) increases when a s tar ascends along the 
AGB lo iofsson et al.' '2002*; 'Wint ers et all 120031: 
iBergeat fc C hcvallier 2005). In this work, we take 
v{oo) calculated by Eq. ([3]) as the standard termi- 
nal velocity of stellar wind. Case 1 is the standard 
model in this work. The results of SSs' population 
are shown in Table [2l 
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Table 1: Parameters of the models of the population of SSs. The first column gives the model number. Columns 
2, 3, 4 and 5 show the mass-loss rate Af , TDU efficiency A, the minimum core mass for TDU and the inter-shell 
abundance, respectively. The detailed descriptions on the above parameters are given in Appendix A. The 6th 
7th and 8th columns show the algorithm of the common envelop e, the terminal velocity and the critical ignition 
mass, respectively. N04 in column 8 means iNelson et al.l ( 2004 



Cases 


M 




A 






Inter-Shell Abundances 


Common 


Envelope 


ti(oo) 


AA/WD 

cnt 


case 1 


Eq. 


(IMt 


Eq.( 


SI 


Eq.(E3l) 


Marigo & Girardi 


r2007'l 


O^ceAcc — 


0.5 


Eq.© 


Eq.gD 


case 2 


Eq. 


(|All| 


Eq.d 


'M 


Eq.(E3l) 


Marigo & Girardi 


f2007') 


O^ccAcc — 


0.5 


Eq.© 


Eq.© 


case 3 


Eq. 


(IMt 


0.5 




Eq.(E3l) 


Marigo & Girardi 


f2007'l 


Q^ccAcc — 


0.5 


Eq.© 


Eq.© 


case 4 


Eq. 


(IMt 


0.75 




Eq.(E3l) 


Marigo & Girardi 


r2no7'i 


Q^ceAcc — 


0.5 


Eq.© 


Eq.© 


case 5 


Eq. 


(lA8l) 


Eq.( 


Ml 


O.58M0 


Marigo & Girardi 


(2007) 


Q^ccAcc — 


0.5 


Eq.© 


Eq.© 


case 6 


Eq. 




Eq.( 


21 


Eq.(R3l) 


Izzard et al. f2004 




Q^ceAce — 


0.5 


Eq.© 


Eq.© 


case 7 


Eq. 




Eq.( 


Ml 


Eq.(R3l) 


Marigo & Girardi 


(2007) 


7 = 1.75 




Eq.© 


Eq.© 


case 8 


Eq. 




Eq.( 


M 


Eq.(R3l) 


Marigo & Girardi 


(2007) 




0.5 


2 ^csc 


Eq.© 


case 9 


Eq. 




Eq.d 


M 


Eq.(R3l) 


Marigo & Girardi 


(2007) 


Q^ccAcc — 


0.5 


Eq.© 


Eq.(Al) of NG4 



Table 2: Different models of the SSs' population. The first column gives the model number according to Table 
[TJ Column 2 shows the birthrate of SSs in the Galaxy. From columns 3 to 5, the number of 0-rich SSs, C-rich 
SSs and all SSs are given, respectively. In the 3rd and 4th columns, the numbers in parentheses mean the 
numbers of 0-rich and C-rich SSs in which we neglect the effects of the mass transfer in binary systems on the 
chemical abundances on their surfaces, respectively. Columns 6 and 7 show the ratios of the number of SSs in 
the cooling phase and the stable hydrogen burning phase to their total number, respectively. Columns 8, 9 and 
10 give the occurrence rates of SyNe (symbiotic novae) with the accreting CO WDs, ONe WDs and total rates. 
The number of SyNe with accreting CO and ONe WDs are shown in columns 11 and 12, respectively. 



Cases 


Birthrate 


Number of SSs ^tcooi 


^stable 


Occurrence 


rate of SyNe (yr 


"^) Number of SyNe 


of SSs (yr^i 


) 0-rich C-rich Total A^totai 


-^total 


COONe 


Total 


CO 


ONe 


1 


2 


3 4 5 6 


7 


8 9 


10 


11 


12 


case 1 


0.119 


4800 (4860) 870 (810) 5670 0.57 


0.25 


3.8 0.2 


4.0 


1000 


4 


case 2 


0.118 


5850 (6070)1240(1020) 7100 0.42 


0.45 


3.4 0.2 


3.6 


930 


3 


case 3 


0.119 


5340 (5390) 220 (170) 5560 0.59 


0.24 


3.7 0.2 


3.9 


940 


3 


case 4 


0.119 


4930 (4970) 840 (800) 5770 0.58 


0.24 


3.9 0.2 


4.1 


1000 


4 


case 5 


0.111 


4330 (4410)1220(1140)5550 0.56 


0.26 


3.8 0.2 


4.0 


1000 


4 


case 6 


0.119 


4740 (4830 930 (840) 5660 0.57 


0.25 


3.8 0.2 


4.0 


1000 


4 


case 7 


0.160 


7900 (7930)1440(1410)9340 0.49 


0.33 


5.3 0.4 


5.7 


1700 


12 


case 8 


0.072 


2370 (2450) 700 (620) 3070 0.28 


0.60 


1.1 0.1 


1.2 


370 


4 


case 9 


0.115 


7830 (7900)1790(1720)9620 0.75 


0.15 


9.8 0.4 


10.2 


1000 


4 
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3.1. Galactic Birthrate and the Number of 
SSs 




Fig. 1. — Initial-final mass relations. Thin and 
thick solid lines are the core masses at the first 
thermal pulse in iHurlev et al. 1 (I2OOOI) and in our 
work (See Eq. (lAip ). respectively. Thin and thick 



dashed lines are the final masses in .Hurley et al 



(|2000l) and case 1 in this work, respectively. The 
dotted lines (from top to bottom) r epresent the re - 
lations from Girardi et al. I (|200d) . iHerwid (|l995l) 
and Weidemann ( 200Cl[ ). respectivelv. The dot- 



dashed line represent the relation from lHan et al.l 
(|1994 ). 



First, we discuss the gross properties of the 
modeled population of SSs as those in Paper I. 
As Table [2] shows, the Galactic birthrate of SSs 
may range from 0.072 (case 8) to 0.160 (case 7) 
yr~^ and their number is from about 3070 (case 
8) to 9620 (case 9). The occurrence rate of the 
symbiotic novae is between 1.2 (case 8) and 10.2 
(case 9) yr~^ and the number of the symbiotic 
novae is from about 370 (case 8) to 1700 (case 7). 
The contribution of the symbiotic novae with ONe 
WD accretors to the occurrence rate is negligible 
in all cases, which is greatly different from that in 
Paper I. The main reason is the different stellar 
evolutions adopted during the TP-AGB phases. 

All assumptions in cases 1 and 8 are respec- 
tively the same as those in case 7 and the stan- 
dard model of Paper I except the stellar evolutions 
during the TP-AGB phase. The greatest differ- 
ence of the TP-AGB phase between this work and 
Paper I is the initial-final mass relations. It is 
well known that the relation is still very uncer- 
tain. Fig. [T] shows the initial-final mass relations 
in different literatur e s. Us ing the mass-loss rate of 



Vassiliadis &: Wood (|l993l) and assum ing the clas- 



sical core mass luminosity relation, IWeidemann 
((2000) gave an initial-final mass relation which is 
plotted in F ig. [H The relation f or the more mas- 
sive stars in iGirardi et all ()2000l ) is steep because 
TDU and hot bottom burning were not included. 
By applying new observational data (including 
NGC2168, NGC2516, NGC2287 and NGC3532) 
and improved theory (stellar evolution calcula- 
tions including new opacities and main sequence 
overshoot; the extensive mass grid of WD cooling 
sequences includi ng mode l s with thick and thin 
He- and H- layer), iHerwid (|l995h gave an initial- 
fi nal mass rel ation which has a less slope than that 
in lGirardi et a l. (2000) for the more massive stars. 
Based on new stellar evolution calculations (in- 
cluding TDU and hot bottom burning) and new 
observational data, Weid emann (2000) showed an 
initial-final mass relation whic h has lower final 
mass than that of Herwig ( 1995h in the more mas- 
sive stars. Using a cri terion for envelop e ejection 
in AGB or FGB stars. iHan et al.l (|l994l) obtained 
an initial-final mass relation which is very sim- 
ilar to that in Weidemann ( 200 0,). P a per I ac- 
cepted the description of IHurlev et al. ( 2000[ ) on 
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TP-AGB evolution. The thin sohd hne and the 
dashed hne in Fig. [T]represent the core mass at the 
first thermal pul se Mc.itp and th e initial-final mass 
relations in [Hurley et al. ( 20001 ) . respectively. In 
this paper, we obtain the final mass by the syn- 
thesis TP-AGB evolutions (Detailed descriptions 
can be see n in Appendix A ). The Mr.Mr, (See 



Eq. (8) in Karakas. Lattanzio fc Polsll2002l ) and 



the final mass are shown by the thick solid and 
dashed lines in Fig. [H re spectively. The rela- 
tion of Hurlev et al.1 (|200Ci ) is steeper than ours 
for Minitiai > 4Mq (especially Minitiai > QMq) 
and gives a higher WD's mass than ours. In Pa- 
per I, we assumed that the ONe WD originates 
from a star with initi al mass being be tween 6.1 



Mq and about 8 M^ (|Pols et al.Ml998D . In this 
work, we still keep the above assumption. The 
masses of ONe WDs in this work are much lower 
than those in Paper I. A small mass of the accret- 
ing WD is unfavorable for producing the symbiotic 
phenomenon. Compared with case 7 and the stan- 
dard model in Paper I, the birthrate, the number 
of SSs and the symbiotic novae decrease in cases 1 
and 8. Especially, the occurrence rate of the sym- 
biotic novae with accreting ONe WDs in case 1 is 
~ 3% of that in case 7 of Paper I, and in case 8 it 
is ~ 5% of that in the standard model of Paper I. 
Therefore, the initial-final mass relation has great 
effects on the SSs' population and especially on 
the symbiotic novae. 

Usi ng the observational d ata in the last 25 
years, llben fc Tutukwl (|l996l ) estimated that the 
occurrence rate of symbiotic novae is ~ 1 yr~^ 
within a factor of 3. As Table. [2] shows, the 
occurrence rate in case 1 is 4 and close to 



the oc currence rate estimated by llben fc Tutukov 
( 19961 ). Under the same conditions, case 7 in Pa- 
per I gave a higher occurrence rate of 11.9 yr~^. 
Obviou sly, an steep initial- final mass relation like 
that in iHiirlev et all (l2000[ ) overestimates the oc- 
currence rate of symbiotic novae. Thi s indicat es 



that the initial-final mass relations fromlHan et al 



()l994l ). IWeidemannl (|2000h and the present paper 
are more reliable than steep ones. However, due 
to the complicity of SSs, the above point of view 
needs supports from further theoretical study and 
more observational evidence. 



3.2. The Effects of Input Parameters 

In this section, we discuss the effects of the 
input physical parameters on the SSs' population 
and the ratio of the number of C-rich SSs in which 
the cool companions are C-rich giants to that of 
0-rich SSs in which the cool companions are O- 
rich giants. In Paper I, we analyzed the effects 
of the physical parameters ( the common enve- 
lope algorithm, the terminal velocity of the stellar 
wind t;(oo) and the critical ignition mass of the 
novae AM^-^ ) on the SSs' population. They 
have a great effect on the SSs' population in our 
simulations, which is similar to those in Paper I. 
Table [5] shows, the common envelope algorithm 
and AA/j^f have weak effects on the ratio of the 
number of O-rich to that of C-rich SSs while the 
v{oo) in case 8 takes an uncertainty within a fac- 
tor of about 1.6. With stars ascending along the 
AGE, the v{oo) in case 8 decreases while more 
C-rich giants are formed. A small t;(oo) is favor 
for the formation of SSs. Therefore, the ratio of 
O-rich to C-rich SSs in case 8 decreases. The ef- 
fects of other physical paramet ers are following. 
Mass- loss rate: Compared with Vassiliadis fc WoodI 
1993 ) (see Eq. (jA8|) ). the mass loss rate of Blockeil 



19951 ) fsee Eq. (jAllI) ) greatly depends on stellar 



luminosity. When the cool giants ascend along the 
AGB, the luminosities gradually increase. There- 
fore, case 2 gives a higher mass-loss rate in the 
later phase during the AGB, which can increase 
the SSs' number of the stable hydrogen burn- 
ing while decrease the SSs' number of the ther- 
monuclear runaways. In general, C-rich giants 
are formed at a later phase of the AGB by the 
third dredge- up. The mass-loss rate in [Blocker 
(|l995l )fSee Eq. (jAlip i enhances the number of 
SSs and the number ratio of C-rich to O-rich SSs 
but decreases the occurrence rate of the symbiotic 
novae and their number. However, the effect is 
weak. 

TDU efficiency A : Comparing cases 1, 3 and 4, 
we find that A has a great effect on the ratio of 
O-rich to C-rich SSs' population. The effect of A 
on stellar evolution has mainly two aspects, i.e. 
the core mass evolution and the changes of the 
chemical abundances in stellar envelope. The av- 
erage dredge-up efficiency in case 1 is higher than 
that in case 3 (A = 0.5) and approximated to that 
in case 4 (A = 0.75). In cases 1, 3 and 4, the SSs' 
populations are almost the same. However, the 
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number ratio of 0-rich to C-rich SSs in case 3 is 4 
times more than the ratio in cases 1 or 4. 
The minimum core mass for T DU M™'": For the 
same i nitial masses, the M™'" of lKarakas. Lattanzio k, 
(|2002l )rSee Eq. jM])) is higher than O.SSMg. In 



rich giants may have high terminal velocity of the 
stellar wind ^(oo ) . Wh en we use the fit of v{oo) 



Winters et al. (l2003ll. we should note that the 



case 1, a star with an initial mass higher than ^ 
2.0Mq can undergo the third dredge-up, but only 
1.5M0 in case 5. Therefore, C-rich SSs are more 
easily produced in case 5. Its effect is within a 
factor about 1.5. 

The inter-shell abundances: As cases 1 and 6 in 
Table [2] show, the two different inter-shell abun- 
dances have i io effe ct on the SSs' populations. In 
Izzard et al. ( 20041 ) . ^^O abun dance of the inter- 
shell r egion is lower than that in lMarigo fc Girardi 
(|2007l ). The number of C-rich SSs in case 6 in- 
creases. However, the effect of the inter-shell 
abundances on ratio of 0-rich to C-rich SSs is 
very weak. 

In short, the mass-loss rate and the inter-shell 
abundances have very weak effect on the SSs' pop- 
ulation and the ratio of 0-rich to C-rich SSs. The 
TDU efficiency A has a very weak effect on the 
SSs' population, while the number ratio of 0-rich 
to C-rich SSs have a great dependence on A. 

3.3. C-rich SSs 



Belczvhski et al. ( 2000[ ) presented a catalogue 
of SSs which includes 188 SSs as well as 28 ob- 
jects suspected of being SSs. According to the 
spectral types of the cool components listed by 
the catalogue, one can find out 5 C-rich SSs in 
176 Galactic SSs and the number ratio of 0-rich 
to C-rich SSs is about 35. In our simulations, 
from cases 1 to 9, the number ratios of 0-rich 
to C-rich SSs are about 5.5, 4.9, 24.1, 5.8, 3.6, 
5.2, 5.6, 3.4 and 4.3, respectively. Our results are 
lower than the above observational ratio. In all 
cases, the ratio of case 3 with a low A=0.5 is th e 



closes t to the ratio observed in iBelczvhski et al 
(I2OOOI) . However , according to the typic a l TP 



19931: 



AGB s ynthesis (IGro enewege n fc de Jond 

Karaka s. Lattanzio k, P ols 200^ Izzard et al I2OO4I : 
Marigo fc Girardill2n07h . A=0.5 may be too low 



for the solar-like stars. In the solar neighborhood, 
the number r atio of M- type giant s to C-type gi- 



ants is ~ 5 in 



Groenewegen (|2002h . which is close 
to our results except for case 3. We consider that 
the small fraction of the observed C-rich giants 
in SSs may result from two respects: (i)The C- 



Polal ample of Winters et al. (|2003h has 65 sources 
but only two of them are C-rich giants. Eq. ([3]) 
may be unsuitable for the C-rich giants. Com- 
paring cases 1 and 8, one can find that v{oo) has 
a great effect on the number ratio of 0-rich to 
C-rich SSs. We may have underestimated v(oo) 
of C-rich giants. (ii)The cool giants in SSs have 
high mass-loss rat es. In cases 1 and 2 we us e the 
mass-loss rate s of IVassiliadis fc WoodI (Il993h and 
iBlocked (Il995l) . respectively. The average mass- 
loss rate of the former is higher than that of the 
later during TP-AGB when the TDU occurs. The 
higher the mass-loss rate is, the more quickly the 
envelope mass decreases. According to Eq. (|A2[) 
in the present paper, the interpulse period Tip in- 
creases with the envelope mass decreasing. A long 
Tip reduces the TDU progressive number and the 
TDU efficiency A. The large mass-loss rate is un- 
favorable for the formation of the carbon stars, 
which can be seen by comparing the number ratio 
of 0-rich to C-rich SSs in case 1 with that in case 
2. iMikolaiewskal (|2007l ) suggested that both the 
symbiotic giants and Miras have higher mass-loss 
rates than single giants or field Miras, respectively. 
Therefore, the predicted number ratio of 0-rich to 
C-rich SSs lower than the observational value may 
result from our underestimating the v{oo) of C- 
rich giants and the mass-loss rates of the cool 
giants in SSs. 

Based on the formation channels, C-rich stars 
are classed into two types. The intrinsic C-rich 
stars mean that the carbon observed in the at- 
mosphere of the TP-AGB stars results from the 
third drudge-up. The extrinsic C-rich stars are 
the giants polluted by carbon-rich matter from 
the former TP-AGB companion. In this work, 
we consider the effects of the mass transfer in bi- 
nary systems on the chemical abundances on their 
stellar surfaces (See iJ2.2p . In order to check the 
above effects, we also carry out the simulations in 
which the mass transfer between the two compo- 
nents can not vary the abundances on their sur- 
faces. The results are shown in the parentheses 
of columns 3 and 4 of Table. [D The ratios of 
SSs with the extrinsic C-rich giants to the total 
C-rich SSs from case 1 to case 9 are 6.8%, 17.7%, 
22.7%, 4.7%, 6.6%, 9.7%, 2.08%, 11.4% and 3.9%, 
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respectively. The ratio is sensitive to the efficiency 
of TDU A, the common envelope algorithm and 
the mass-loss rate. According to the discussions 
in WS.21 we know that a high mass-loss rate occurs 
in case 2 when carbon abundance is increased by 
the TDU. The extrinsic C-rich giants in case 2 are 
formed more easy than those in case 1 . A low A in 
case 3 can prevent the formation of C-rich stars in- 
cluding intrinsic and extrinsic C-rich stars. The 
high ratio in case 3 results from a small number of 
C-rich stars. Having undergone the a-algorithm 
common envelope evolution, the orbital period of 
binary system becomes to be several percent of 
that before common envelope evolution. With 7- 
algorithm in case 7, post-common-envelope sys- 
tems are wider than with the a-algorithm in case 
1 and this facilitates a symbiotic phenomenon, al- 
lowing more stars to evolve further along the FGB 
or AGB before the second RLOF (Paper 1) . In gen- 
eral, post-common-envelope systems hardly form 
the extrinsic C-rich stars but can form intrinsic 
C-rich stars. Therefore, the ratio of SSs with the 
extrinsic C-rich giants to the total C-rich SSs in 
case 7 is very low. In most cases, the fraction of 
extrinsic C-rich SSs in total C-rich SSs is lower 
than 10%. This indicates that the pollution of the 
carbon rich materials from the former TP-AGB 
companions on the cool components is weak so 
that most o f C-rich giants in SSs are intrinsic C- 
rich stars. Mikolajewskal ( 2007 ) showed that all 
C-rich giants observed in SSs are intrinsic C-rich 
stars, which is in agreement with our results. 



from the observations. ISchild. Bovle &: Schmid 



3.4. 



VS. [C/H] of the Cool Giants 

in SSs 



(11992) used the synthetic spectr a calculated for 
M giants bv iLazaro et al. ( 1991 ) in which they 
analyzed M giants an d obtained a low mean 
[C/H]=-0.64 ± 0.29. iLazaro et al.l (Il99lh sug- 



gested that the sta ndard mixing named as the 
first dredge-up (See llben fc Renzinil 119831 ) is in- 
sufficient to explain atmospheric abundances in 
M giants. In our work, we use the standard mix- 
ing which only includes the first dredge-up. ^^C 
abundance on stellar surface is reduced by approx- 
imatel y 30% after the first dredge-up during the 
F GB |lben fc Renzinilll983f ). 

Gratton et all (|200d ) determined Li, C, N, O, 



Na and Fe abundances and ^^C/^'^C isotopic ratios 
of 62 metal-poor field giants. They suggested that 
there are two distinct mixing stages along the red 
giant branch evolution of small mass field stars: (i) 
The first dredge-up, that is, the standard mixing; 
(ii) The second mixing episod e which is also called 
the th ermohalinc mixing in Charbonnel fc ZahnI 
(I2OO7I) when the giants have a brighter luminosity. 
The thermohalinc mixing can decrease ^^C abun- 



dances. Charbonnel fc Zahni ( 2007h calculated 



several evolution models of a O.9M0 star and ob- 
tained the abundances wh ich are consistent wit h 
the observational data in ICratton et all (|2000l ). 
Recently, lEggleton etld] (1200^ |2007t) confirmed a 
possible mechanism for such non-canonical mix- 
ing. Using 3D-modeling of a low-mass star at 
the tip of red giant branch, they found that the 
molecular inversion can lead to an efficient mixing 
and Pop I stars between 0.8 and 2.0 Mq develop 
i^ C/"C ratios of 14.5±1.5 The a verage of [C/H] 



Schild. Bovle fc Schmid I (|l992l ) is about -0.54 



Using th e infrared spectra of SSs,| Schild, Bovle fc Schmiff| d the average ^"^C/^^C is about 16. It range 



(|1992| ) and ISchmidt fc Mikolajewsk a (2003) gave 
the C abundance and ^^G/^^G isotopic ratio of 
the cool giants in 13 SSs with ±0.3 errors. Us- 
ing high-resolution near-infrared spectra and the 
method of the standard local thermal equil ibrium 
analys is and atmospheric models, iSchmidt et al.l 
( 20061 ) calculated the abundance of the symbi- 
otic star CH Cyg with the errors less than 0.3 
dex for all the elements. Their observational data 
are shown in panels (10) of Fig. ^ where [C/H] 
means the relative abundances to the solar, that 
is, [C/H] = logC/H-logC0/H0. The left panel 
of Fig. [2] shows that ^^C/^^G and [C/H] of the 
cool giants in SSs are much higher than those 



in our simulations is between about -0.05 ( case 
3) and 0.15 (case 4) for [C/H] and between about 
23.4 (case 3) and 37.2 (case 5) for ^^G/^^G, respec- 
tively. The main reason of the above disagreement 
is that the thermohalinc mixing is not included in 
our simulation. In order to show the effects of 
the thermohalinc mixing, we artificially decrease 
^^G abundance of all giants with initial masses 
< 2Mq by a factor of 2.5 after the first dredge- 
up. The results are given in the right panel of 
Fig. [21 The averages of [C/H] are between -0.31 
(case 3) and -0.102 (case 5) and the averages of 
i2(-;y'i3(^ are between 12.3 (case 3) and 19.5 (case 
5). Our results are then close to the observations 
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1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 
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1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 
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Fig. 2. — Gray-scale maps of the chemical abundance ratios of log ^^C/^'^C versus [C/H]. The gradations of 
gray-scale correspond to the regions where the number density of systems is, respectively, within 1 - 1/2, 
1/2 - 1/4, 1/4 - 1/8, 1/8 - of the maximum of gipg lac/^^caiog [c/h] ' ^^'^ blank regions do not contain 
any stars. Solid lines mean log ^^C/^'^C=log [C/H]. The left panel of Fig. [2] shows the chemical abundance 
ratios of the cool giants in SSs. The right panel of Fig. [2] gives those of the cool giants in SSs, where ^^C 
abundance of all giants with initial masses < 2Mq is artificially decreased by a factor of 2.5 after the first 
dredge-up. In the pan el of observations, the filled big -stars, th e empty squares and trian gle represent the 
o bservational val u es in Schild. Bovle fc Schmid ( 1992 ). those in lSchmidt &: Mikolaiewskal |2003( ) and those 
in Schmidt et al. ( 20061 ). respectivelv. 



13 



of lSchild. Bovle fc Schmidi (|l992h . 

We believe that the thermohaline mixing can 
occur and give significantly different resuhs al- 
though the observational samples of ^^C/^'^C vs. 
[C/H] of the cool giants in SSs is small. However, 
due to poor knowledge on the thermohaline mix- 
ing, we can not carry out a detailed simulation. In 
the following sections, we do not discuss its effects 
any more. 

3.5. Chemical Abundance of Symbiotic 
Nebulae 

In this subsection, we discuss the chemical 
abundances of the symbiotic novae. In order to 
clarify the nature of the symbiotic nebulae, we 
need the observational abundances of novae and 
planetary nebulae . In this paper, the no va abun- 
dances come from Livio fc Truran ( 19941 ) in which 
they analyzed the abundances of 18 classical no- 
vae with a factor of ~ 2-4 uncertainty. The 
planetary nebulae abundan ces ar e obtained from 
|Pottasdi_& Bernard-Salas ( 20061 ) in which they 
gave the chemical abundance of 26 planetary neb- 
ulae within a 30% uncertainty in the abundance 
determination. 

3.5.1. C/N vs. 0/N 

Based on UV data, iNussbaumer et aL I (|l988l) 



presented C/N and 0/N abundance ratios of 24 
symbiotic nebulae. The errors in the logarith- 
mic (basis of 10) abundance ratios remain within 
0.18. In their sample, there are at least three 
symbiotic novae V1016 Cyg, HM Sge and HBV 
475. V1016 Cyg and HBV 475 had their outbursts 
around 1965 while HM Sge brigh tened by 5 mag 
in 1975. 1 Schmid fc Schild I (|l990h calculated their 
abundance ratios within an error of 30%. Their 
C /N an d O /N abundance r a tios a re different from 
those in INussbaumer et al.l (|l988r ) within a factor 
of 30%. According to lUE observations, HM Sge 
is still evolving towards higher excitation whereas 
V1016 Cyg and HBV 475 seem to have reached 
their maximum excitation (jNussbaumer fc Vogel 
1990|). The outburst of PU Vul began late in 



1977. IVogel fc Nussbaumeil (|l992[ ) gave the C/N 
and 0/N abundances ratios in the early nebular 
phase. The observational results and our simu- 
lated results are shown in Fig. [3l 

The left and right panels (1)— (9) of Fig. dgive 



the distribution of the relative CNO abundances 
of the cool companions and the ejected materials 
from the hot components during symbiotic novae, 
respectively. Panels (10) in Fig. [3] show the ob- 
servational ratios of C/N vs. 0/N of symbiotic 
nebulae, novae and planetary nebulae. 

The left panel of Fig. [3] shows that we can cut 
the abundance ratios of C/N vs. 0/N of the cool 
companions in SSs into three regions, i) The re- 
gion of C /N > O /N represents that the cool com- 
panions have undergone the deep TDU and have 
not undergone the hot bottom burning. The ini- 
tial masses of the cool companions should be be- 
tween about 2.0 or 1.5(case 5) Mq and 4.0 Mq. 
Their abundances are similar to those of C-stars 
and planetary nebulae; ii) The top region of C/N 
< 0/N represents that the cool companions have 
undergone inefficient TDU or have not undergone 
TDU. These stars have initial masses lower than 
2.0 or 1.5 (case 5) Mq. Their abundances are sim- 
ilar to those of M-stars; iii) The bottom region of 
C/N < 0/N represents that the cool companions 
undergo TDU and the hot bottom burning which 
turns to "N. Their initial masses are higher 
than 4.OM0. 

In the right panel of Fig. [3l the distribution of 
C/N vs. 0/N of the ejected materials from the hot 
companions during symbiotic novae is composed 
of two regions. The top region represents symbi- 
otic novae with accreting ONe WDs. The bottom 
region represents those with accreting CO WDs. 
Compared wi th the observations, ou r results agree 
with those in iLivio fc TruTa^ (|l994[ ). 

From the left panel of Fig. [3l we find that the 
relative CNO abundances of the cool companions 
in our simulations are close to those of the obser- 
vational symbiotic nebulae. However, comparing 
more carefully, one can find that many observed 
SSs lie in the transition from normal M giants with 
initial mass lower than AMq (no hot bottom burn- 
ing) to the super giants with higher initia l mass. 
This agrees with INussbaumer et al. I (Il988l) . If we 
compare Fig. [Sjwith the observations, we find that 
the abundance ratios of most of observational sym- 
biotic nebulae are between those of the cool giants 
and the ejected materials from the hot compan- 
ions during symbiotic novae. This suggests that 
the compositions of the symbiotic nebulae may 
be modified by the ejected materials from the hot 
companions. 
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Fig. 3. — Gray-scale maps of the chemical abundance ratios of log C/N versus log 0/N. The gradations of 
gray-scale correspond to the regions where the number density of systems is, respectively, within 1 - 1/2, 
1/2 - 1/4, 1/4 - 1/8, 1/8-0 of the maximum of giogc/N3iogO/N ' ^^'^ blank regions do not contain any stars. 
Solid lines mean log C/N=log 0/N. The left and right panels of Fig. [3] show the chemical abundances of the 
giants during the symbiotic stage and the ejected materials from the hot companions during the symbiotic 
novae, respectively. In the panels of observations, the filled b ig-stars, the filled triangl es and the squares 
re present the observational values of the symbiotic nebulae in Nussbaumer et al. ( 19881 ) . those of PU Vul 
Vogel fc Nu s sba umer (1992) (including Civ and Njy or excluding them) and those of Cyg, HM Sge and 



HBV 475 in fSchmid fc Schild I (|l99Cl( ). respectively. The circles in the left and right pan e ls (10 ) of Fig. 
m repres ent the observat i onal v alues of the planetary nebulae in Pottasch fc Bernard-Salad ( 20061 ) and the 
novae in lLivio fc Trurarj (|l994f ). respectively. 
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As the panels (10) of Fig. [3] show, the abun- 
dance ratios of the symbiotic novae Pul Vul and 
HM Sge are hke those of the novae while the ratios 
of V1016 Cyg and HBV 475 are similar to those of 
the giants. Pul Vul and HM Sge are still staying 
on the early phase of the nova outbursts whereas 
V1016 Cyg and HBV 475 may be staying o n the 
decline phase. Vogel fc Nussbaumer ( 1992h pre- 
dicted that the chemical abundances of Pul Vul 
would become closer to those of red giants in the 
future. Therefore, we believe that the composi- 
tions of the symbiotic nebulae have been modi- 
fied by the ejected materials from the hot com- 
panions during symbiotic novae. Can the modifi- 
cation only last ton like our assumption in ^2AV 
According to Table [2l most of SSs are in the cool- 
ing phase of the symbiotic novae and the stable 
hydrogen burning phase. The symbiotic novae are 
only 18.3% (case 7) and 10% (case 9) of total SSs. 
However, most of the observed symbiotic nebulae 
have mixing CNO abundances, which means that 
the above modification should be common in all 
SSs. We showed the C/N and 0/N abundance ra- 
tios of several typical SSs which stay in different 
phases in order to support the above the view of 
point. 

Declining phase: AG Pegasi is currently in 
decline from a symbiotic nova outburst which be- 
gan in about 1850. Based on multi-shell radio 



emissi on from AG Pegasi in lKennv. Tavlor fc Seaquist 



(|l991l ). the nebula can be divided into three parts. 
The outer nebula is thought to be the remnant 
of the 1850 outburst. The intermediate nebula 
may represent the pre-eruption mass loss from 
the cool component. Using a new colliding winds 
models includin g the effects of orbital motion, 
Kennv &: Tavloii (|2007l) suggested that the inner 
nebula originates from the colliding winds and has 
an expansion velocity of about 50 km s~^. It takes 
about 1000 yr for the inner nebula to reach the re- 
gion of the outer nebula. In general, the timescale 
of the cooling phases tcooi is several hundred years. 
Therefore, during the cooling phase after the ther- 
monuclear runaways, the ejected materials from 
the hot components have an effect on the chemi- 
cal compositions of the symbiotic nebulae. 

Multiple outbursts: ISokoloski erall (|2006l ) 



shell burning triggered by a disk instability. Dur- 
ing the burst occurring in 2000-2002, a shell of 
material blown from the surface of the WD is 
found. The blown material may result from an 
optically thick wind. Both Z And and CH Cyg 
are SSs which have undergone multiple outbursts. 
Since 196 4, CH Cyg has displa yed a number of 
outbursts. Schmidt et al. ( 20061 ) showed both the 
photospheric features of the cool giant and the 
nebular emission lines in CH Cyg. The average 
C/N and 0/N ratios of the cool giants in CH 
Cyg are respectively 1.6 and 4.0, which are typical 
abundance ratios of M giants. The average ratios 
of the symbiotic nebulae in CH Cyg are respec- 
tively 0.57 and 2.0, which are similar to those of 
super giant stars and are between those of giants 
and novae. The significant discrepancy of C/N 
and 0/N ratios derived from these two compo- 
nents suggest that the chemical abundances of the 
nebulae in CH Cyg are modified by the ejected 
materials from the hot components. 

Quiescent phase: Sy Mus lacks any outburst 
activity and it is in the sta ble hydrogen burning 
phase ( Mikolajewskal l2007l ). iNussbaumer et al 



(|1988[ ) showed that the C/N and 0/N ratios in 
its nebulae are 0.97 and 1.9, respectively. These 
ratios are similar to those of super giant stars and 
are between those of giants and novae. This means 
that the chemical abundances of the nebulae of Sy 
Mus in the quiescent phase are modified by the 
ejected materials from the hot components. 

Therefore, we suggest that it is quite common in 
SSs that the chemical abundances of the symbiotic 
nebulae are modified by the ejected materials from 
the hot WD surfaces. 

3.5.2. Ne/0 vs. N/0 

According to optical data. lLuna fc Costal (|2005h 
calculated He abundance and the relative abun- 
dances N/0, Ne/0 and Ar/0 of 43 symbiotic neb- 
ulae. We first discuss Ne/0 vs. N/0 . The mean 
errors in the logarit hmic (basis of 10) abundance 
ratios are ± 0.1. IPe Freitas Pacheco fc Costal 



discussed the outburst of symbiotic system Z And 
occurring in 2000-2002. They believed that the 
outburst in 2000-2002 results from the enhanced 



(11992) showed the ratios of N/0 and Ne/0 in 5 
symbiotic nebulae within an error of 0.15. Using 
the observational data, we note that the observa- 
tional results for same SSs in different literature 
are quite different. The ratio s of N/0 and Ne/0 
of V 1016 Cyg in Dc Freit as Pacheco fc Costa 
(jl992l ) differ by a factor of 4 from those in 
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Fig. 4. — Gray-scale maps of the chemical abundance ratios of log Ne/O versus log N/0. Solid lines mean 
log Ne/0=log N/0. The left and right panels of Fig. |4]show the chemical abundances of the giants during 
the symbiotic stage and the ejected materials from the hot companions during the symbiotic novae, respec- 
tively. In the panels of the observations, the filled bi g-stars, the filled tria ngles and the squares represent 
th e observational values of the symb iotic nebulae in iLuna &: Costal (|2005f ) , th ose of 5 symbiotic nebula e 
in iDe Freitas Pacheco fc Costal (|l992l ) and those of V1016 Cyg and HM Sge in I Schmid fc Schild I (|l990l ). 
respectively. The cir cles in the left and right panels ( 10) of Fig. [H repr esent the observationa l values of the 
planetary nebulae in Pottasch fc Bernard-Salad ( 20061) and the novae in Livio fc Truran (1994), respectively. 
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Schmid fc Schild I (|l990ll . 

Fig. m gives the distributions of the abundance 
ratios of Ne/0 vs. N/0 in all cases and observa- 
tions. The left and right panels of Fig. 2] show 
Ne/0 vs. N/0 of the cool giants in SSs and the 
ejected materials from the hot components during 
symbiotic novae, respectively. The right panel of 
Fig. Ulshows that our results of the ejected materi- 
als from the hot companions during the symbiotic 
novae are less scattered than the observations of 



the novae in iLivio fc TruranI (|1994| ). There are 
mainly two reason s. One reason is too smal l nu- 
merical samples in Kqvetz fc Prialni^ ( 1997 ) and 
Jose fc Hernanj ( 1998h . Another reason is a high 
uncertainty in observations with a factor of 2-4. In 
the right panel of Fig. [4l the distributions of Ne/O 
and N/0 are obviously cut into two regions. The 
left-up region of higher Nc/0 originates from the 
ejected materials from the hot companions dur- 
ing the symbiotic novae with accreting ONe WDs 
and the middle-down region of lower Ne/0 corre- 
sponds to the ejected materials from the hot com- 
panions during the symbiotic novae with accreting 
CO WDs. 

Panels (10) of Fig. 0] show that the abun- 
dance ratios of Ne/0 vs. N/0 in symbiotic neb- 
ulae are basically between those in the planetary 
nebulae and the novae. Our simulations also sup- 
port this result. Therefore, the distributions of 
Ne/0 vs. N/0 reconfirm that the symbiotic neb- 
ulae are modified by the ejected materials from 
the hot components. 

3.5.3. He/H vs. N/0 

Luna fc Costal (|2005l ) showed the distributions 



of the observational log 0/N vs. He/H for 43 
symbiotic nebulae. The mean error of log 0/N 
is about 0.1 and the mean error of He/H is w ithin 
~ 0.03. ICosta fc de Freita s Pachecol (|l994l ) gave 
those of 5 symbiotic nebulae within an error of 
0.038. 

As the panels (10) of Fig. [5] show, there are 
many symbiotic nebulae whose helium abundances 
are much higher than those of the observational 
planet ary nebulae in Pottasch fc Bernard-Salad 
( 20061 ). The enhancement of the helium abun- 



stellar surface of the cool giants in all SSs are 
basically lower than 0.2 in our simulations. The 
enhancement of the helium abundance in the sym- 
biotic nebulae should originate from the ejected 
materials during the symbiotic novae. The ob- 
servational Hc/H ratios of the classical novae in 
fLivio fc TruranI (jl99j) range from 0.07 to 0.40. 
The observational distribution of He/H vs. log 
N/0 also support that the symbiotic nebulae are 
modified by the ejected materials from the hot 
components. Unfortunately, as the right panels 
(l)-(9) of Fig. El show, our simulations do not 
offer high He/H ratios for the ejected materials 
from the hot components. A possible explanation 
is: if the mass-accretion rates of WD accretors 
are higher than the critical value Mwd (See Eq. 
(24) in Paper I), the nova outbursts are weak and 
a helium layer may be left on the surface of WD 
accretors. We call these novae as weak symbiotic 
novae (Paper I). If the mass- accretion rates of WD 
accretors are lower than Mwd , the nova outbursts 
are so strong that the WDs are eroded. Therefore, 
no helium layer is left. These novae are called as 
strong symbiotic novae. A typical symbiotic star 
can undergo dozens of nova outbursts (Paper I). 
For the weak symbiotic novae, a helium layer may 
be left on the surface of WD accretors after each 
outburst. When the next outburst occurs, the 
remnant he lium can be dredged up and ejected 
( Iben et al ][l992). \ n addition, most WDs have 
been traditionally found to belong to one of the 
following categories: those with a hydrogen-rich 
atmosphere (the DAs) and those with a helium- 
rich atmosphere (the DBs). DA white dwarfs 
constitute about the 80% of all observed WDs. 
They have a thin hydrogen layers whose mass 
is lower than about 10~'*Mwd- Also, there is 
a helium layer (their mass is about 10~^Mwn\ 



unde r the hydrogen layer (Althaus fc Benvenutol 



19981 ). Ilben fc Tutukovl (|l985h calculated the dis 



tributions of helium, carbon and oxygen through 
the surfaces layers of a 1.05 M© DB WD com- 
ing from a case B mass-transfer event. They 
found the amount of helium near the surface is 
only a bout 10^^ Mq. Gonzalez Perez fc Metcal"^ 



dance in the stellar surface mainly depends on the 
third dredge-up. However, the left panels (1)- 
(9) of Fig. [5] show that the He/H ratios in the 



([20071) tested the evolution of the DB white dwarf 
CD 358. They suggested that binary evolution 
describes better CD 358 and obtained the helium 
layer with a mass lO^^'^^AfyvD- Therefore, when 
the first nova outbursts occurs, the hydrogen shell 
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Fig. 5. — Gray-scale maps of the chemical abundance ratios of He/H versus log N/0. The left and right 
panels of Fig. [5] show the chemical abundances of the giants during symbiotic stage and the ejected materials 
from the hot companions during the symbiotic novae, respectively. In the panels of the observations, the 
filled big-stars , the s quares and the triangles represe nt the observational values o f the s ymbiotic nebulae in 



iLuna fc Costa ( 2005), those of 5 symbio t ic neb ulae in lCosta &: de Freitas Pacheco ( 1994 ) and those of V1016 



Cyg and HM Sge in Schmid fc Schild ( 1990l ). respectively. The circ les in the left and right panels (1 0) of 
Fig. repre sent the observa t ional values of the planetary nebulae in Pottasch fc Bernard-Salad ( 20061 ) and 
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flash will penetrate inward to the region where the 
abundance of H is ~ 0.01 and a certain am ount of 
the helium is dredged up ( Iben et al]ll992h . In the 
above situations, helium is overabundant. Our cal- 
culations may underestimate helium abundances 
of the ejected materials from the hot components 
during the symbiotic novae. 

The helium o f the symbiotic nov a e V10 16 Cyg 



and HM Sge in I Schmid fc Schild I (|l99d) is not 
overabundant. They may have undergone several 
strong symbiotic novae so that the WD accretors 
have no the helium layer. Their helium should be 
similar to our simulated result. 

We suggest that the helium enhancement of the 
symbiotic nebulae stems from the modification of 
the ejected materials from the WDs with a helium 
layer. In our simulations, the ratios of the SSs be- 
ing in the weak symbiotic novae and their decline 
phases to total SSs are between 25% (case 2) and 
38% (case 9). The helium element of these sym- 
biotic nebulae should be overabundant. The ra- 
tio of the symbiotic nebulae whose He/H is higher 
than 0.2 to the tota l of the symbiotic nebulae in 
Lima fc Costal \200^ is 33.3%, which are in agree- 
ment with our results. 

4. Conclusion 

We have performed a detailed study of the 
chemical abundances in SSs and drew several im- 
portant conclusions. 

1. The initial- final mass relation has great ef- 
fects on SSs' population and especially on the 
symbiotic novae with massive WD accretors. A 
steep initial-final mass relation result in an over- 
estimated occurrence rate of symbiotic novae. 

2. The number ratios of 0-rich SSs to C-rich 
SSs in our simulations are between 3.4 and 24.1, 
and they are sensitive to TDU efficiency A and the 
terminal velocity of stellar wind v(oo). Our simu- 
lation may have underestimated the terminal ve- 
locity of stellar wind ^(oo) of C-rich giants and the 
mass-loss rate of the cool giants in SSs. The num- 
ber ratio of SSs with extrinsic C-rich cool giants 
to all SSs with C-rich cool giants is between 2.1% 
and 22.7%. The TDU efhciency A, the common 
envelope algorithm and the mass-loss rate have a 
great effect on it. 

3. Comparing our ^^C/^^C vs. [C/H] of the 
cool giants in SSs with those of the observations, 



we infer that thermohaline mixing in low-mass 
stars should exist. Its effect on the chemical abun- 
dances is very significant. 

4. The distributions of 0/N vs. C/N, Ne/0 vs. 
N/0 and He/H vs. N/0 of the symbiotic nebulae 
indicate that it is quite common in SSs that the 
nebular chemical abundances are modified by the 
ejected materials from the hot components. 

5. Helium abundances in the symbiotic nebu- 
lae during the symbiotic novae are determined by 
whether the WD accretors have a helium layer or 
not in their surfaces. If they have a helium layer, 
helium is overabundant. If the WD accretors have 
only undergone strong symbiotic nova outbursts, 
they have no helium layer. The helium in the sym- 
biotic nebulae ( like V1016 Cyg and HM Sge ) are 
not overabundant. 
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A. Synthesis TP-AGB Evolution 



Stellar evoluti on from the z e ro ag e main sequence up to the first thermal pulse is dealt with in the rapid 
evolution code of iHurlev et alj ( 20001 ). The changes of the chemical abundances on the stellar surface during 
the giant branch phase (the first dredge up) and ea rly asymptot ic gian t bra nch ph a se (th e second dredge 
up) can be represented by simple fitting formulae in Izzard et al. ( 20041 ) and llzzard (2004). After the first 
thermal pulse, we use a synthetic model for TP-AGB. 



A.l. The initial abundances 



The initial abundances {i.e. of zero age main sequence stars) are taken from lAnders fc: Grevesse (Il989h 
for Z = 0.02. The following shows the initial abundances by mass fractions: 



iH=0.68720; 
i2C=2. 92293x10-3 

"N=8.97864xl0"'^: 
1^0=8.15085x10-3 
20Ne=2. 29390x10-3; 



4Hc=0.29280 ; 

i3c=4.10800xl0-5 

i5N=4.14000xlO"^ 

1^0=3.87600x10-6 

22Ne=1.45200xlO- 



A. 2. Core mass at the first thermal pulse 

Using the rapid evolution code of Hurley et al. ( 2000l ). we can obtain the stell ar mass (Mitp) at the first 
thermal pulse. The core mass at the first thermal pulse, Mc.itp is taken from iKarakas. Lattanzio fc Fold 
(|2002l) : 

MciTP = [-piiMo - P2) + P3]f + {PiMo + P5)il - /), (Al) 

/ Mq —PQ \ 

where / = (1 + )~^, Mn is the initial mass in solar unit and the coefficients Pi, P2j Pa, Pi, Ps, Pe ^nd 

p^ are in Table 6 of iKarakas. Lattanzio &: Fold (|2002h . 

A. 3. Luminosity, radius, interpulse period, and evolution of core mass 



We use the prescrip t ions o f llzzard et al.l (j2004l ). The luminosity is taken as the value calculated from 
Eq.(29) in llzzard et all (|2004f) . We define the radius J? as L = AwaR'^T^ff, where a is the Stefan-Boltzmann 
constant a nd T^s is th e effec tive temperature of the star. The radius is taken the value calculated from 
Eq.(35) in llzzard et al.l (|2004l ). The interpulse period Tip is: 



logio(T-ip/yr) = a28(Me/M0 - bsg) - 10^'^* 



10^ 



0.15A^ 



(A2) 



where the third dredged-up (TDU) efficiency A is defined in ^A.Al and the coefficients are: 

= -3.821, 
&28 = 1.89 26, 

C28 = -2.080 - 0.353Z + 0.200(Afc„v/MQ + a- 1.5), 
d28 = -0.626 - 7O.3O(Afc,iTp/M0 - C)(AAfc/Af0), 

where Afcnv represents the envelope mass, a is the mixing length parameter and equals to 1.75, C = 
log(Z/0.02), and AAfc is the change in core mass defined as AAfc = Mc — Mc^np. 

A. 4. The minimum core mass for TDU and the TDU efficiency 

TDU can occur only f or sta rs above a certain core mass Af™'". iGroenewegen fc de Jond ( 1993[ ) took 
Af™'" as a constant 0.58. IKarakas. Lattanzio fc Polsi (,2002) found that A^T™™ depends on stellar mass and 
metallicity and gave a fitting formula by 



ai + a2Mo + asM^ + a^M;^, 



(A3) 
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where coefScients ai, 02, 03 and 04 are shown in Table 7 of lKarakas. Lattanzio k. Pold ()2002h and Mq is the 
initial mass in sol a r unit . According to the observed carbon luminosi ty function in the Magellanic clouds, 
Marigo fc Girardil (|2007t ) considered that the M™'" predicted by Kar akas. Lattanzio fc Pols (,2002) is high. 



In this work, we carry out yarious num erical simulations (See Table [T]). 

(i) Like'Groenewegen & de Jong" figgi), M^=O.58Af0; 

(ii) Following Karakas, Lattanzio & Vo\^ (|2002l ). we take Eq. ((X3|l as M^'"'"; 

It should be recalled th at M™™=Mc_itp if stellar initial mass Afinitiai > 4Mq ( Karakas. Lattanzio fc Pol^ 
I2OO2I : llzzard et al.ll2004h or M™ < Mc4TP- 

The TDU efficiency is defined by A = '^am"'' ' where AMdrcd is the mass brought up to the stellar surface 
during a thermal pulse. A is a very uncertain parameter. Karakas. Lattanzio fc Poli (2002) showed a relation 
of A: 

X{N) = A,nax[l - exp(-iV/A^,)], (A4) 
where A gradually increases towards an asymptot ic Amax with N (th e progressive number of thermal pul 



sation) incre asing. N-^ is taken from Eg . (49) in Izzard et al.l ( 20041 ) which re produces the results for iVr 



in Tab le 5 of lKarakas. Lattanzio fc Pold ( 2002 ). Amax is given by (See Eq.(6) in Karakas. Lattanzio fc Poli 
(|2002l) 1: 

bi + + 63^3 ^^^^ 



1 + 64M3 

where coefficients are shown in Table 8 of lKarakas. Lattanzio fc Pold ( 2002 ) . Groenewegen fc de Jone* (Il993l ) 
took A as a constant 0.75. In order to test the influences of A on our results, we carry out numerical simulations 
with yarious A (See Table HI . 

CDFoUowing 'Karakas. Latt anzio fc P"ols ('2002^, A=Eq. (|A4l) : 

(ii) Following Gro enewegen fc de Jon g (1993), A=0.75; 

(iii) Simulating a small TDU efficiency, A=0.5. 

A. 5. Inter-shell abundances 

During every thermal pul se, the dredged mass AMdred mixes into stellar envelope . Based on the nu- 
cleosynthesis calculations by iBoothrovd fc SackmannI (|l988r ). iMarigo fc Girardil (|2007r ) gave the fit of the 
abundances of ''lie, ^^C and ^^O in the inter-shell region: 
For AMc <0.025Mq: 

inter^(-4jjg) ^ o.95 + 400(AMc)2 - 20.0AMc 
irAer xii-iQ) = 0.03 - 352(AMc)2 + 17.6AMc 

inter js^|l60) = -32(AMc)2 + 1.6AMc. 

For AMc >0.025Mo: 

mter^(-4jjg) ^ o.70 -I" 0.65(AMc - 0.025) 

inter js^(12Qj ^ q 35 _ 0.65(AAfc " 0.025) 

intcr^^ieg) ^ 0.02 - 0.065(AAfc - 0.025). 



However, there is some debate of the exact composition in the inter-shell region, llzzard et al.l (|2004l ) gave 
some fits of inter-shell abundance, in whic h inter-shell elements include ^He, ^^C, ^^O and ^^Ne. Details can 
be seen from jj3.3. 3 in Izzard et al. I (l2004h Their models did not obtain high values of inter-shell ^^O such 
as 2% reported bv lBoothroyd fc SackmarmI (il988,). The typical inter-shell abundances (5 M©, Z = 0.02) are 
4He=0.74, i2c=0.23, 1*^0=0.005 and 22Ne=0.02. 

In our work, we carry out simulations with various inter-shell abundances. 



A. 6. The hot bottom burning 

If the hydrogen envelope of an AGB star is sufficiently massive, the hydrogen burning shell can extend 
into the bottom of the convective region. This process is called as hot bottom burning(IIBB). For HBB, 
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we us e a treatment similar to that in Groenewegen fc de Jond (1993). For the model of Iben fc Renzini 
(1983), Groenewegen fc de Jond (1993) gave the most suitable parameter for the fraction of newly dredged 
up matter exposed to the high temperatures at the bottom of envelope /hbb — 0.94, the fraction of envelope 
matter mixed down to the bottom of the envelope /b ur — 3 X 10 ^, and the exposure time of matter in the 
region of H BB tHBB = 0.0014Tip. The temperature at the bottom of convective envelope Tbco is given by 
Eq. (37) of llzzard et all (|2004). 

For the drudged up mass, the amounts of material added to the envelope are: 



A^He = 
A12C = 

A13C = 
A14n = 
Ai6o = 



inter 



X4AAfdr 



[(1 - /hBb)'"*"''Xi2 



Jo 



/hbb 

^HBB 

r /hbb 

^HBB 

[/hbb rtHBB 



vHBB 
^12 



L ^hbb 



Jo 
Jo 



X 



hbb 

13 
vRBB 
^14 



{t)dt]AMdred 
{t)dt]AMdr-ed 



(A6) 



[(1 - /hBb)'"*'='-Xi6 



Jo ^r^Wdi]AMdred, 

where X^^^{t) are calculated by the way of C la yton's CNQ bicyc le('Clavton"l983'). All details of Cla yton's 
CNO bicycle can be seen from Clavton ( 1983[ ): Groenewegen fc d c Jong (1993); I zzard et al. ( 2004h . The 
CNO bicycle can be splited into CN cycle and ON cycle. The timescales in ON cycle are many thousands of 
times longer than those required to bring the CN cycle into equilibrium. Even in the most massive AGB stars 
undergoing vigorous HBB, the ON cycle never approaches equilibrium. Therefore, the effects of HBB mainly 
turn ^^C into ^^N and the abundance of ^^O is not changed much. In c alculating Clayton's CNO bicycle, the 
density at the base of the convective envelope is given by Eq. (42) of Izzard et al.' ('200 4|) and the analyti c 
expressions for the nuclear reaction rates in Clayton's CNO bicycle are taken from Caughl an fc Fowlerl (jl988 ). 
The initial conditions of X^^^{t) are X^^^{t = 0) ='"*"' X for ^^C and ^^O, while X™^(t = 0) = for 
"C and "N. 

After every thermal pulse, the chemical abundances of stellar envelope X 



are 



X°'^Afenv(l - /bur) + AX 



^HBB 



Mr 



AMdrcd 



(A7) 



old 



j^ncw 

where AX is given by Eq. (|A6p and the initial conditions of X{t) are X{t — 0) — X 
A.7. Mass loss 

The mass-loss rate of the cool giant during AGB phase has a great effect on the population of SSs and 
the chemical evolution of the stellar surface. We conside r two l aws of the mass-loss rate in this work. 
(i)A mass-loss relation suggested by Vassiliadis fc Wood ( 1993t ) based on the observations, given as 



logM = -11.4 + 0.0123(P- 100 max(Af/Af0 - 2.5,0.0)), 
where P is the Mira pulsation period in days given by 

logP = -2.07+ 1.94 log(i?/i?0) - 0.901og(M/Mo). 
When P > 500 days, the steady super-wind phase is modeled by the law 

M(MQyr"^) = 2.06 x 10"^^^^^^^^, 

where Woo is the terminal velocity of the super-wind in km s~^. We use Woo = 15 km in this pape r. 

(ii)Based on t he simulations of shock-driven winds in the atmospheres of Mira-like stars in iBediinI (1988), 
Blockeii (jl995fl gave a mass-loss rate similar to Reimers' formula: 



(A8) 
(A9) 

(AlO) 



M 



4.83 X 10-^M-2.i^2.7^^^^^._^^^^^^ 



(All) 
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where MRcimcrs is given by Eq.fJTl) but 77 = 0.02 (jStanchffe fc Jeffervll2007l ). 

On other stcUar evolutionary phase, the mass-loss rates are given by Reimers' formula ( Reimer j 1975 ): 



T ¥t 

M = -4.0 X 10"^^?? Moyr"\ 

M 

where L, i?, M are the stellar luminosity, radius and mass in solar units and we use 77=0.5. 



(A12) 



A. 8. Comparison with previously published yields by other models 

In order to test our synthetic model of asy mptotic giant branch stars, we compa r e the stellar y ields 



in our mode l with previously published yields of Ivan den Hoek fc Groenewegei 
llzzard et al. (2004). According to the definition of stellar yield in llzzard et al.l 
yield of isotope j by 




Marigo (20011) and 



we write the total 



(A13) 



where r(Mi) is the entire lifetime of a star with initial mass Mi, is the current mass-loss rate, Xj and 
Xj refer to the current and initial surface abundance of the isotope j, respectively. 

Figure |6| shows the stellar yields of Ivan den Hoek &: Groenewegei] (Il997l) [for t h e models of Z = 0.02, 
jyAGB = 4.0 and tohbb = 0.8 (See Table 17 of Ivan den Hoek fc Groenewegenlll997l ) 1. iMarigol (|200l[ )(for the 
models of Z = 0.019 and t he mixing length parameter a = 1.68), llzzard et al.l ( 20041 ) [for the models of 
Z = 0.02 (See Table D2 of llzzard et alll2004t )] and ours. In our models, all physical parameters are the 
same as those in case 1 except that we change binary systems in to single star . As F igure El shows, the stellar 
yields of our models in genera l are i n agr eement with t hose in llzzard et al.l (j20041) or lie between those of 
van den Hoek fc Groenewegenl ([19971) and llzzard et al.l (j2004[ ) although our models produce more ^^N, and 



less -^''O at massive stars. 
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Fi g. 6. — Total stellar yields vs. the init i al st ellar masses. T riangles, circl e s and squares represent the models 
of van den Hoek fc Groenewegen ( 1997 ). lMarigo. (,200lh and Izzard et al. ( 20041 ). respectively. Pluses are our 
results. 
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